Time-dependent density-functional approach for exciton binding energies 



V. Turkowski 1 - 2 and C. A. Ullrich 1 

'Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211 
2 Department of Physics and NanoScience Technology Center, 
University of Central Florida, Orlando, FL 32816 
(Dated: February 28, 2009) 

Optical processes in insulators and semiconductors, including excitonic effects, can be described 
in principle exactly using time-dependent density- functional theory (TDDFT). Starting from a lin- 
earization of the TDDFT semiconductor Bloch equations in a two-band model, we derive a simple 
formalism for calculating excitonic binding energies. This formalism leads to a generalization of 
the standard Wannier equation for excitons, featuring a nonlocal effective electron-hole interaction 
determined by long-range and dynamical exchange-correlation (XC) effects. We calculate excitonic 
binding energies in several direct-gap semiconductors, using exchange-only and model XC kernels. 
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The elementary model of Wannier excitons in insula- 
tors views them as bound electron-hole pairs which sat- 
isfy a hydrogen-like Schrodinger equation 
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Here, m r is the reduced electron-hole effective mass, 
e is the electron charge, e is the static dielectric con- 
stant of the material, and <f> and E are the excitonic 
wave functions and binding energies (from now on we 
set U = e = 1). Eq. ([1]), also known as Wannier equation 
0, produces a Rydberg series of discrete energy states 
below the conduction band edge, and a redistribution of 
oscillator strength in the optical spectrum around the 
band edge which is qualitatively described by the Elliott 
formula |3|. Excitonic effects are important for a large 
variety of optical processes in organic and inorganic ma- 
terials and nanoscale systems 

Eq. ([!} can be derived from the semiconductor Bloch 
equations within the time-dependent Hartree-Fock ap- 
proximation, using a dielectrically screened Coulomb in- 
teraction [3, 6]. It is well known that time-dependent 
Hartree-Fock with bare Coulomb interaction leads to 
very poor optical spectra of materials, with strongly over- 
bound excitons. A more rigorous ab-initio treatment of 
excitation processes in insulators and semiconductors, in- 
cluding correlation-induced screening, can be developed 
using many-body Green's function techniques such as the 
GW/Bethe-Salpeter equation 0. 

Time-dependent density-functional theory (TDDFT) 
8] has recently emerged as an alternative, computation- 
ally convenient approach to electronic excitation pro- 
cesses in materials 

IS El E3- In linear -response 
TDDFT, excitation energies can be calculated in prin- 
ciple exactly [H], EH , provided the exchange-correlation 
(XC) kernel / xc (r,r',o;) is known. In Refs. 0, Q, an ap- 
proximate / xc was constructed from many-body Green's 
functions, whereas Ref. [Io| uses an exact-exchange 
(EXX) approach, including a cutoff in wavevector space 



which mimics screening of the Coulomb interaction 14 1 . 
These studies have established that TDDFT is capable of 
accurately describing excitonic effects in solids, although 
one needs XC functionals that go beyond the more com- 
mon ones such as the adiabatic local-density approxima- 
tion (ALDA) Q. The resulting agreement with experi- 



mental data is excellent [ll|, but the technical effort is 



not significantly less than for standard many-body ap- 
proaches. 

The purpose of this paper is to develop a formally 
much simpler TDDFT treatment of excitonic effects in 
solids. Rather than calculating complete optical spec- 
tra, our goal is more modest, namely a method that 
directly yields excitonic binding energies, similar to the 
Wannier equation (fTJ) . Starting from a TDDFT version 
of the semiconductor Bloch equations [l5| . we derive an 
effective electron-hole interaction which explicitly shows 
how long-range XC effects are essential for exciton forma- 
tion. Our simplified treatment not only provides phys- 
ical insight into the way excitonic effects are treated in 
TDDFT, but also provides a straightforward way of test- 
ing approximate XC functionals. 

Time- dependent Kohn-Sham formalism for solids. - In 
TDDFT, the electron dynamics of a solid is described by 
the time-dependent Kohn-Sham (KS) orbitals ^>jk(r,t), 
where k is the wave vector and j is the valence band index 
(only the time evolution of the initially occupied states 
is considered). The system is assumed to start from the 
ground state, 4' J k(r,to) = ^jktr). The KS Bloch func- 
tions and band structure follow from 
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where V\ a t is the crystal lattice potential (within the 
Born-Oppcnhcimer approximation) and Vjj and V~ c are 
the static Hartree and XC potentials. 

Since the V'jk(r) form a complete set for each k, we 
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can expand the time-dependent KS orbitals as follows: 

^• k (r,t)=^4'(f)^k(r) ) (3) 
l 

where the summation runs over all valence and conduc- 
tion bands, including continuum states. Eq. ((3J is ap- 
propriate if we assume the system to interact with an 
electromagnetic field in dipole approximation. We define 
the density matrix p l j£(t) — (^(^[c^ (t)]* , whose equa- 
tion of motion is 
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where we dropped the time-dependent external field 
term, since the excitations we are interested in can be 
viewed as eigenmodes of the system. Here, 8V^ and 
5V™k denote the linearized dynamical Hartree and XC 
potentials. In a periodic insulating solid, the Hartree 
term only gives rise to the so-called local field correc- 
tions, which do not affect excitonic binding We will 
therefore only keep the XC contribution in the following. 

Fourier transformation of Eq. §9§ and the correspond- 
ing equation for p^(t) leads to 
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with initial condition Pj£(to) — §ji8 m i. The matrix ele- 
ments of the TDDFT Hamiltonian are 

= e^5 lm + E(t)dL m + + V&(t) , (5) 

where f2 is the volume of the lattice unit cell, E(£) is the 
electric field amplitude, and d^ 71 are the dipole matrix el- 
ements. Vn(t) = Vn(t) — Vjj denotes the dynamic part of 
the Hartree potential, and similar for XC. Self-consistent 
solution of Eq. (j4]) , with the time-dependent density 
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where ef is the Fermi energy, is equivalent to solving the 
time-dependent KS equations for the solid, and is thus 
in principle exact. 

Two-band model and excitons. - To study optical exci- 
tation processes near the band gap, a two-band model is a 
reasonable and widely used approximation. We consider 
one valence and one conduction band, v and c, assumed 
to be nondegenerate (see [lfj} for a discussion of band de- 
generacy). The index j of the density matrix p l J£(t) refers 
to v and will be dropped in the following. Eq. (j4|) yields 
the TDDFT semiconductor Bloch equations for the two 
independent components p k and p k c 03 : 
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Notice that p£ v + p k c = 1 and p k c = p%>*. In Ref. [l5(, 
Eqs. ([7]) and © were evaluated in the time domain for 
ultrafast pulsed excitations. Here, we are interested in 
excitonic binding energies, and we linearize Eq. (|8j): 



}pl c (t)-6V^(t)~6V^ k (t), (9) 
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d 3 r' < k (r)^- k (r)/ xc (r,r',w) 
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X Cq( r ') f /'nq 

and the q-summation runs over the first Brillouin zone. 
Eqs. (|10p and (fTTj) can be cast into an eigenvalue prob- 
lem for the excitation energies ui. Since / xc is in general 
frequency-dependent, the eigenvalue problem is nonlin- 
ear. The solutions are the exact excitonic binding ener- 
gies within the two-band model. 

Let us carry out a further simplification. Since typi- 
cal excitonic binding energies are much smaller than the 
band gap, i.e., uj + lu^ ui — tu^, we can ignore the pole 
at negative oj (which is equivalent to the Tamm-Dancoff 
approximation and boldly set p k c = 0. This leads to 
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Eq. ([13]) is the equivalent for extended systems of the 
well-known single-pole approximation of linear-response 
TDDFT [12]. For finite atomic or molecular systems, 
the single-pole approximation only involves two discrete 
levels. Here, it involves two entire bands, which clearly 
shows the collective nature of excitonic effects. 

We point out that Eq. (fT3|) yields excitonic binding 
energies relative to the conduction band edge, which can 
be accurate even if the band gap itself is not. 

TDDFT Wannier equation. - Our next goal is to de- 
rive a real-space equation for the excitonic binding ener- 
gies. p k u is a periodic function in reciprocal space, with 



Fourier transform p(R, to) = Ek e p^fo), where R 
is a direct lattice vector. Similarly, we define 

V eh (K,R',u>) = ^ e - lk ' R F k c ™ c ( w ) e ^ R '. 
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From the point of view of a Wannier exciton, which ex- 
tends over many lattice constants, R can be approxi- 
mated as a continuous variable. We assume a direct band 
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gap material, and use approximate parabolic dispersions 
with conduction and valence band effective masses m e 
and m v , and reduced electron- hole effective mass m^ 1 = 
to" 1 + m^ 1 . This yields the TDDFT version of the Wan- 
nier equation |T]), 

— Eb,% / d 3 r'V eh (r,r',u)) Pi (r')=0, 

space 

(15) 

featuring a nonlocal, frequency-dependent electron-hole 
interaction V e h(r, r', w), where lo = E^ s + E^i, and £* s 
is the KS band gap. The ith excitonic binding energy 
Ej, t i is measured with respect to the KS conduction band 
edge, and the pi (r) are the analog of the excitonic wave 
functions </>(r) of Eq. (JTJ. 

XC kernels. - The effective electron-hole interaction in 
TDDFT, and thus the excitonic binding energies, depend 
crucially on the approximate XC kernel. In the following, 
we shall implement several simple frequency-independent 
XC kernels and test their performance in our formalism. 

The exchange-only ALDA kernel is given by 



/ x ALDA (r,r') = -[9^(r)]- 1 /3 (5(r _ r ' ) 
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where no(r) is the equilibrium electron density. / ALDA 
belongs to the class of ultra-short-range kernels; the sim- 
plest of them is 
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where A is a positive constant. Such kernels have been 
used with some success in contact exciton models |lll|. 

An approximation of exact-exchange TDDFT [l(| , the 
Slater exchange kernel, is given by [121 ] 
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|r - r'|rao(r)no(r') 
This kernel exhibits some degree of long-range behavior 



17 1 , but not the ultra-nonlocality (~ 1 /q 2 in momentum 



space) of the exact / xc @, EH • This long-range contri- 
bution (LRC) can be explicitly taken into account using 
the following model kernel [181 ]: 
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where a is again an adjustable parameter. 

Results and Discussion. - We have tested our TDDFT 
approach for excitonic binding energies, Eq. (fT3"|) . for 
the zincblende materials GaAs and /3-GaN and for the 
wurtzite materials a-GaN, CdS, and CdSe. The Bloch 
functions for the conduction and heavy-hole valence 
bands were obtained from LDA band structures calcu- 
lated with the plane-wave pseudopotential code ABINIT 
[l9j |. We used 512 k-points in the first Brillouin zone 
for all materials. Out of these, there are 10 independent 
points for GaAs and /3-GaN and 20 for a-GaN, CdS, and 



TABLE I: E^ p and £° latcr : lowest direct excitonic binding 
energies (in meV) for selected III-V and II- VI compounds, 
from experiment H2 and from Eqs. [fI5]). 1(12)1 with / x ilater . 
The parameters A and a/4n (in a.u.) are fitted to reproduce 
E^ p using Eqs. JT2J with f^° ntact and /^ C RC . 





A 


Q/47T 


rpSlatcr 




GaAs 


1.68 


0.12 


17.4 


3.27 


/3-GaN 


4.23 


0.55 


28.7 


26.0 


a-GaN 


2.03 


0.91 


11.8 


20.4 


CdS 


6.28 


1.83 


7.9 


28.0 


CdSe 


4.84 


1.19 


8.3 


15.0 



CdSe, which determines the dimension of the eigenvalue 
problem (TT3"|) . Recent Bethe-Salpeter calculations of ex- 
citonic binding energies used much higher k-point densi- 
ties close to the zone center 2(| 21 1; we performed con- 
vergence checks of our k-point sampling rates and found 
them to be sufficiently accurate for our simple model. 

As expected, the ALDA does not produce any bound 
excitons. Results for the other three XC kernels and ex- 
perimental binding energies of the lowest direct excitons 
are presented in Table U The contact and LRC kernels, 
([17]) and (fl~9j) , contain adjustable parameters which can 
be tuned to reproduce the experimental exciton binding 
energies. The required values of the parameters A and a 
are found to be of similar order as in Ref. [l8( . 

The contact and LRC kernels only yield a single exci- 
tonic bound state [l8| . This is generally the case for static 
XC kernels that are local in reciprocal space, i.e., have 
the form / xc (q)- The Slater XC kernel (fT5)) does have 
some degree of nonlocality in reciprocal space, but we 
found that it only produces a single excitonic state, like 
the local kernels. To obtain an excitonic Rydberg series 
one needs an XC kernel that has a sufficiently strongly 
nonlocal form or is frequency-dependent [23|, l24j . 

Looking at the results obtained with /x later , we find 
excitons that are overbound by 14 meV in GaAs and by 
2.7 meV in /3-GaN. This overbinding is what one would 
expect from an unscreened exchange-only approach (elec- 
tronic screening can be viewed as a correlation effect). On 
the other hand, /Slater approaches a constant for q — > in 
homogeneous systems [171 ] . whereas the full EXX / x be- 
haves as 1/q 2 [l0( . This would suggest that / x later has a 
somewhat weaker effective electron-hole interaction than 
full EXX. This trend seems confirmed in the wurtzite 
materials whose calculated excitonic binding energies are 
significantly below experiment. 

Additional insight is provided by comparing the 
electron-hole interaction V e t for the different XC kernels 
under study. Fig. [1] shows F e /i(r, 0) for GaAs along the x 
direction (due to the finite sampling in k-space, V e h can 
only be reliably calculated within the range of about one 
unit cell). In ALDA, the interaction is close to zero and 
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in wurtzite. There are theoretical arguments in favor of 
both trends, which suggests a need for more systematic 
studies of the Slater exchange kernel in solids. 

In conclusion, our simple approach for excitonic bind- 
ing energies is a promising method to test XC kernels in 
solids. It can be extended in a straightforward way to 
deal with spin-dependent excitations (triplet excitons), 
with more sophisticated XC kernels, or to include more 
bands. Work along these lines is in progress. 

This work was supported by NSF Grant DMR- 
0553485. We thank Angel Rubio, Lucia Reining, and 
Claudia Ambrosch-Draxl for useful discussions. 



FIG. 1: Effective electron-hole interaction V e h{v, 0) for GaAs 
and different XC kernels, plotted along the x direction, where 
a is the lattice constant. The parameters A and a for the 
contact and the LRC XC kernels are given in Table [I] 



thus too shallow to lead to any excitonic binding. The 
other XC kernels produce stronger electron-hole interac- 
tions, where for GaAs the contact and LRC models are 
less attractive than the Slater approximation. 

Conclusion. - We have presented a simple method 
to calculate excitonic binding energies using TDDFT. 
The main idea, restricting the dynamics to the high- 
est valence and the lowest conduction band, is similar to 
the single-pole approximation for excitation energies [12j . 
Our derivation was based on the TDDFT semiconduc- 
tor Bloch equation; an alternative starting point could 
be the Casida formalism of linear-response TDDFT [l]| , 
formulated for periodic systems [2oT |. The resulting sim- 
ple eigenvalue equation in momentum space, Eq. (|13j) . is 
readily diagonalized to yield the excitonic binding ener- 
gies. Transformation into real space leads to the TDDFT 
analog of the Wannier equation for excitons, and shows 
that the effective electron-hole interaction is nonlocal. 

The quality of the results depends crucially on the ap- 
proximation used for / xc (r, r', lu). It is well known that 
local and semilocal approximations, such as the ALDA, 
do not produce any excitons. There exist sophisticated 
parameter-free XC kernels Q that are capable of repro- 
ducing experimental optical absorption spectra very ac- 
curately, including bound excitons [24j |. but with sub- 
stantial computational cost. 

If only particular aspects of the optical spectrum of 
a material are required such as, for instance, the lowest 
bound exciton, simple static XC kernels can be a con- 
venient alternative. The contact and the LRC kernels 
behave quite similarly in the sense that they produce a 
single excitonic peak. A detailed analysis was given in 
Ref . [231 ] , and we find the same behavior in our two-band 
approach. The parameter-free Slater exchange-only ker- 
nel also produces a single exciton, which was found to 
be overbound in zincblende materials, and underbound 
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